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metal-metalloid  glasses,  mainly  because  they  are  the  oldest  and  best  studied 
group.  Significant  progress  has  also  been  made  in  the  study  of  metal-metal 
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metal-metalloid  glasses. 
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by  structural  imperfections.  As  discussed  in  the  second  section,  defects  in 
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concentration  and  motion  of  individual  defect  sites  that  produce  local  shear. 
Again,  this  implies  a certain  limitation  in  scope,  since  significant  development^ 
concerning  transient  flow,  anelasticity,  internal  friction  and  embrittlement 
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I.  INTRODUCTION 


This  paper,  if  it  is  to  stay  within  reasonable  limits,  cannot 
possibly  be  as  general  as  the  title  indicates.  Therefore,  rather  than 
trying  to  cover  all  the  recent  developments  in  the  study  of  structure  and 
flow,  an  attempt  will  be  made  to  provide  a link  between  the  two  subjects: 
the  emphasis  will  be  on  those  structural  aspects  that  are  important  for 
flow  and  on  a description  of  the  flow  process  that  is  as  consistent  as 
possible  with  the  structural  findings. 

In  a first  section,  recent  developments  in  the  study  of  the  structure 
of  amorphous  alloys  are  reviewed.  The  discussion  is  limited  to  the 
transition  metal-metalloid  glasses,  mainly  because  they  are  the  oldest 
and  best  studied  group.  Significant  progress  has  also  been  made  in  the 
study  of  metal- metal  and  rare  earth-tr ansition  metal  alloys,  but  these 
developments,  although  quite  important,  have  not  been  included.  It  is 
shown  how  the  new  scattering  techniques  have  confirmed  earlier  ideas 
about  the  high  degree  of  short  range  order  in  the  metal- metalloid  glasses. 

In  amorphous  systems  as  well  as  in  crystals,  atomic  transport  is 
governed  by  structural  imperfections.  As  discussed  in  the  second  section, 
defects  in  amorphous  metals  are  probably  more  diffuse  than  in  crystals 
and  can  therefore  best  be  described  as  individual  sites  with  small  per- 
turbations of  the  local  short  range  order. 

In  the  third  section  this  is  applied  to  problems  of  steady  state 


plastic  flow.  First,  the  experiments  are  reviewed  by  means  of  an 
empirical  deformation  mechanism  map.  Then,  an  attempt  is  made  to 
describe  the  basic  processes,  homogeneous  and  inhomogeneous  flow, 
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from  a unified  point  of  view  involving  the  concentration  and  motion  of 
individual  defect  sites  that  produce  local  shear.  Again,  this  implies  a 
certain  limitation  in  scope,  since  significant  developments  concerning 
transient  flow,  anelasticity,  internal  friction  and  embrittlement  have  not 
been  included. 

f 
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II.  THE  STRUCTURE  OF  TRANSITION  META  L- META LLOID 

ALLOY  GLASSES 

1.  Background 

Three  years  ago,  around  the  time  of  the  previous  conference,  the 
status  of  our  knowledge  about  the  structure  of  amorphous  metals  had 
been  thoroughly  reviewed  by  Cargill.  » Up  to  that  time,  information 
about  the  atomic  scale  structure  of  metal-metalloid  alloy  glasses  derived 
to  some  extent  from  density  measurements  and  indirectly  from  observ- 
ations of  crystallization  and  glass  transition,  but  mainly  from  single 
wavelength  X-ray  scattering  experiments,  which  yield  radial  distribution 
functions  that  are  dominated  by  the  metal-metal  pairs  because  of  the 
larger  concentration  and  scattering  factors  of  the  metal  atoms. 

These  studies  had  demonstrated  conclusively  the  inadequacy  of 
microcrystalline  structural  models  and  the  importance  of  models  based 
on  the  dense  random  packing  (DRP)  of  hard  spheres.  While  the  DRP  of 
single  size  spheres  was  clearly  the  appropriate  model  for  pure  amorphous 
metals,  its  suitability  for  describing  the  structure  of  amorphous  alloys 
was  still  a matter  of  debate,  mainly  because  of  the  lack  of  precise  in- 
formation on  the  position  of  the  metalloid. 

3 

Polk  had  suggested  that  the  metal  atoms  form  a DRP  skeleton, 

with  the  metalloids  occupying  the  larger  holes  in  this  structure.  However, 

since  all  these  holes,  as  described  by  Bernal,^  are  actually  too  small  to 

5 

accommodate  the  metalloids,  Polk  later  modified  his  proposal  by 
appealing  to  the  distortion  of  the  idealized  holes  in  the  actual  DRP  and  by 
allowing,  if  necessary,  some  additional  relaxation  of  the  hole  shape, 
resulting  in  a metal  skeleton  that  is  less  densely  packed  than  the  DRP. 
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The  qualitative  implications  of  this  model  are  that  a metalloid  atom 
(i)  is  surrounded  by  only  metal  atoms  and  (ii)  has  a lower  coordination 
number  than  a metal  atom.  The  same  requirements  were  the  basis  of 
the  binary  dense  random  packings  by  Sadoc  et  al.  , ^ which  had,  however, 
limited  usefulness  for  detailed  comparison  because  of  the  small  size  of 
the  clusters . 

2.  Recent  Experimental  Developments 

The  main  experimental  advances  of  the  past  few  years  are 

scattering  techniques  which  yield  information  about  the  position  and 

7 

environment  of  the  metalloid  atoms.  Sadoc  and  co-workers  used  neutron, 
polarized  neutron  and  X-ray  scattering  to  obtain  the  three  partial  radial 
distribution  functions  for  COg^P^.  They  used  neutron  and  X-ray 
scattering  to  obtain  the  Pd- Pd  and  Pd-Si  partial  radial  distribution 
functions  in  ^dg^Sijg  and  to  establish  the  absence  of  Si-Si  nearest  neigh- 

g 

bors.  Waseda  and  co-workers  made  use  of  the  large  anomalous  dis- 
persion terms  in  the  X-ray  scattering  factor  for  wavelengths  near  the 
absorption  edge  of  the  scattering  element.  They  applied  this  in  the  case 
of  Ni  by  using  Cu  Kqt  and  Co  Kor  radiation.  By  combining  these  data  with 
regular  measurements  away  from  the  absorption  edge  (Mo  Ko),  they 

obtained  the  three  partial  radial  distribution  functions  for  Nior.Pofl. 

oU  ZU 

9 

Hayes  et  al.  used  the  extended  X-ray  absorption  fine  structure  (EXAFS) 
of  Ge  to  determine  the  nearest  neighbor  environment  of  this  metalloid  in 
Pd^gGe^^.  The  nearest  neighbor  structural  data  resulting  from  these 
investigations  are  listed  in  Table  1.  For  comparison,  the  same  data  for 
crystalline  intermetallic  compounds  of  similar  compositions  are  also 
listed. 


TABLE  1:  Nearest  neighbor  structural  data  for  amorphous  and  crystalline  transition  metal  (Me)  - metalloid  (X)  alloys. 
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The  most  important  conclusion  of  these  studies  is  the  confirmation 

of  earlier  speculation  that  the  degree  of  chemical  short  range  order  in 

these  systems  is  very  high;  each  metalloid  is  surrounded  by  metal 

atoms  only.  Comparison  with  similar  crystalline  compounds,  shows  that 

the  metalloid  coordination  numbers  and  the  average  metal- metalloid 

distance  are  similar.  Besides  from  analogy  with  crystalline  ordering, 

the  high  degree  of  chemical  ordering  in  the  glassy  alloys  had  been  antici- 

3 5 15 

pated  on  the  basis  of  model  building,  ’ thermal  behavior  and  compo- 
sition dependence  of  some  properties,^  and  recently  an  attempt  has  been 

17 

made  to  explain  it  in  terms  of  pseudo- potential  theory. 

The  metalloid  radius  can  be  estimated  by  subtracting  the  metal 
radius,  which  for  a transition  or  noble  metal  is  probably  close  to  the 
Goldschmidt  radius,  from  the  average  metal-metalloid  nearest  neighbor 
distance.  The  results  are  listed  in  Table  1 and  illustrated  on  Figure  1. 

As  had  been  known  for  the  crystalline  materials,  ^ ^ and  had  been  anti- 
cipated for  the  glasses,  the  observed  metalloid  radii  are  smaller  than 

the  metalloid  Goldschmidt  radii.  The  observed  radius  for  P is  close  to 

its  tetrahedral  covalent  radius,  as  had  been  pointed  out  in  Rundqvist's 
1 8 

review  on  the  crystalline  phosphides.  The  observed  radii  for  Si  and 

Ge  in  both  crystals  and  glasses,  are  appreciably  smaller  than  their 

respective  covalent  radii.  This  difference  has  also  been  noticed  by 
1 8 

Rundqvist,  but  no  explanation  for  it  has  been  offered.  It  must  be  con- 
cluded that  it  is  difficult  to  choose  a priori  an  appropriate  metalloid 
radius  on  chemical  grounds.  There  is  probably  little  significance  to  the 

P- radius  being  close  to  the  tetrahedral  covalent  one,  since  recent  NMR- 

20 

measurements  on  NiP  alloys  have  shown  that  the  binding  is  metallic  in 


nature. 
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Fig.  1:  Metalloid  radius  vs.  Goldschmidt  radius  for  selected  glassy 
(dots)  and  crystalline  (triangles)  alloys.  Data  from  Table  1; 
data  for  crystalline  Co^P  and  Pd3P  from  Ref.  18;  Goldschmidt 
radii  from  Ref.  10;  covalent  radii  from  Ref.  19. 
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There  seems  to  be  disagreement  between  the  experimental  deter- 
minations of  the  spread  in  metal-metalloid  nearest  neighbor  distances. 

0 

Using  the  EXAFS  technique,  Hayes  et  al.  measure  an  R.  M.  S.  full  width 

of  less  than  0.  2 A,  similar  to  what  is  observed  in  crystals.  This  is 

substantially  less  than  what  is  measured  with  the  other  scattering 

2 1 

techniques,  even  if  the  appropriate  corrections  are  applied  to  the 

listed  FWHM.  Until  this  question  is  resolved,  it  remains  difficult  to 

draw  structural  conclusions  from  this  type  of  observations. 

The  average  metal-metal  nearest  neighbor  distances  in  both  classes 

and  crystals  are  slightly  larger  than  the  metal  Goldschmidt  diameter,  as 

shown  in  Table  1 and  on  Figure  2.  (Some  earlier  data  from  composite 

radial  distribution  functions  have  been  added,  since  they  are  dominated 

by  the  metal-metal  pairs.  ) Figure  2 shows  how  the  ratio  of  this  distance 

to  the  metal  Goldschmidt  diameter  increases  with  increasing  metalloid 

1 23 

content.  As  has  been  pointed  out  before,  ’ this  reflects  the  distortion 
caused  by  inserting  the  metalloids  in  the  pure  metal  skeleton.  It  is 
interesting  to  note  that  this  distortion  is  roughly  twice  as  large  for 
crystals  than  for  glasses  of  the  same  composition.  This  higher  degree 
of  distortion  could  possibly  be  the  result  of  the  additional  constraints 
imposed  on  the  metal  skeleton  by  the  crystal  symmetry  requirements. 

In  order  to  have  a density  slightly  higher  than  the  glasses,  the  distance 
increase  in  the  crystals  is  accompanied  by  an  increase  in  coordination 
number. 

24 

Recently,  Turnbull  has  analyzed  the  densities  of  glassy  and 
crystalline  transition  metal- metalloid  alloys  for  a large  variety  of  compo- 
sitions. As  illustrated  by  Figure  3 for  the  phosphides,  he  found  that  the 
metalloid  partial  gram-atomic  volume  V^.  was  equal,  within  ±10% 


METAL- METAL  N 


at  % METALLOID 


2:  Change  in  the  ratio  of  the  metal-metal  nearest  neighbor  distance 
to  the  metal  Goldschmidt  diameter  with  metalloid  content.  The 
data  for  the  pure  metals,  the  series  of  NiP  glasses,  PdgpSi^Q 
and  CoggP22  are  from  Cargill's  review,  ref.  1;  the  series 
of  FeP  glasses  is  from  ref.  22;  the  rest  is  from  Table  1 (dots: 
glasses;  triangles:  crystals). 
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VMe  (cm  3/g~atom) 


Fig.  3.  V : partial  gram-atomic  volume  of  P in  some  glassy  (dots) 

P and  crystal-line  (triangles)  Me-P  alloys.^ 

V : gram-atomic  volume  of  pure  Me  at  300  K.  After 
Me  Turnbull,  ref.  24;  some  data  from  ref.  22  have  been  added. 
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error  limits,  to  the  gram-atomic  volume  of  the  pure  metal  V°,  (which 

Me 

was  taken  to  be  equal  to  the  partial  gram- atomic  volume  of  the 

metal  in  the  alloy).  This  scaling  relation  indicates  that  the  size  of  the 
metalloid  in  these  alloys  is  not  an  intrinsic  property  of  the  element,  but 
is  primarily  determined  by  a characteristic  spacing  in  the  metal  host 
structure.  However,  when  the  observed  metalloid  radius  and  the  metal 
Goldschmidt  radius  are  compared  (see  Fig.  1),  the  data  are  scattered  and 
no  scaling  is  apparent.  It  is  possible  that  a more  appropriate  choice  of 
size  parameters  might  improve  this. 

3.  Recent  Model  Studies 

5 

As  discussed  above,  Polk's  reason  for  modifying  his  original 
proposal  was  that  the  idealized  shapes  of  the  holes  that  make  up 

4 

Bernal's  DRP  are  too  small  to  accommodate  the  metalloids.  At  the 
time,  this  was  based  on  analogy  with  the  crystalline  metalloid  environment. 
It  is  confirmed  by  the  measurements  of  the  metal-metalloid  distance  in  the 
glasses,  as  shown  by  the  'interstitial  size  ratio'  in  Table  1.  This  ratio, 
defined  as  the  metal- metalloid  distance  over  the  metal  Goldschmidt  dia- 
meter, is  no  smaller  than  0.  87,  compared  to  a center-to- vertex  over 
edge  length  ratio  of  0.  82  for  the  largest  idealized  hole.  1 

The  possibility  remained,  however,  that,  since  the  actual  holes  in 
the  DRP  are  distorted  (in  the  sense  that  their  edge  lengths  can  differ  by 
20%),  they  could  accommodate  larger  interstitial  spheres  than  the 
idealized  holes.  In  the  light  of  this  possibility,  F :>st  has  recently 

investigated  the  hole  structure  of  Finney's  mechanically  built  and 
27 

Bennett's  computer- built  DRP.  He  observed  that,  with  spheres  within 
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a fixed  cutoff  distance  being  nearest  neighbors,  these  DRP  structures 

2! 

cannot  be  completely  described  as  stacking  of  the  five  Bernal  polyhedra, 
but  a number  of  new  shapes  with  four-  and  five-edged  faces  have  to  be 
included,  especially  to  describe  the  larger  holes.  While  confirming  the 
predominantly  tetrahedral  character  of  the  DRP  structure,  he  has 

obtained  statistics  for  the  large  holes  which  are  quite  different  from 

4 .45 

Bernal's.  A similar  observation  has  been  made  by  Whittaker. 

25 

For  this  reason,  Frost  has  made  an  independent  calculation  of 
the  exact  number  and  size  of  all  the  interstitials  in  the  DRP,  without 
reference  to  the  hole  structure.  This  was  done  by  considering  the  inter- 
stitial spheres  in  all  sets  of  four  DRP  spheres  (since  the  center  of  the 
largest  interstitial  sphere  is  always  equidistant  between  four  DRP  sphere 
centers)  and  eliminating  the  ones  that  overlap  with  a DRP  sphere  or  a 
larger  interstitial.  The  results  are  shown  in  Figure  4.  The  main  con- 
clusion of  this  study  is  that,  although  the  holes  are  distorted  in  the  DRP, 
they  are  not  appreciably  larger  than  the  idealized  ones:  the  largest 


center  to  center  distance  is  0.  86  and  in  order  to  accommodate  25  inter- 


stitials per  100  DRP  spheres  the  interstitial  site  ratio  must  be  smaller 
than  0.  74.  A study  of  the  interstitials  in  the  Bennett  DRP  led  to  similar 
conclusions.  Fitting  the  metalloids  of  Table  1 in  a DRP  skeleton,  there- 
fore, requires  some  relaxation  of  the  structure.  The  remaining  question 
is  whether  this  relaxation  involves  only  local  adjustments  around  the 
metalloid,  or  whether  it  requires  a reconstruction  of  the  entire  model. 


Some  insight  on  this  question  can  probably  be  obtained  from  the 
recent  development  in  constructing  dense  random  packings  of  binary 
mixtures  of  different  size  spheres.  Using  the  solution  of  the  Percus- Yevick 
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CENTER  TO  CENTER  DISTANCE,  R 

26 

Fig.  4.  Interstitial  sphere  statistics  for  Finney’s  DRP,  compared 
to  Bernal's**’  * original  statistics.  (After  Frost,  ref.  25).  R 
is  the  distance  between  the  center  of  the  interstitial  sphere  and 
the  center  of  one  of  the  four  equidistant  nearest  DRP  spheres. 
Unit  length:  DRP  sphere  diameter. 
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equation,  Weeks  V has  calculated  the  partial  pair  correlation  functions  for  a 

number  of  binary  hard  sphere  mixtures  with  various  compositions  and  sphere 

size  ratios.  This  approach  however  fails  to  take  into  account  the  chemical 

interaction  between  the  two  kinds  of  atoms.  This  makes  it  inappropriate  to 

describe  the  structure  of  metal- metalloid  glasses,  as  demonstrated  by  its 

prediction  of  a large  number  of  metalloid-metalloid  nearest  neighbors. 

The  high  degree  of  chemical  ordering  was  an  explicit  element  of  the 

computer-built  binary  hard  sphere  DRP's  of  Sadoc  et  al.  ^ This  approach  has 

30  31 

recently  been  refined  by  Boudreaux  and  Gregor  ’ who  have  built  much 

larger  models  and  have  paid  special  attention  to  the  problems  of  isotropy 

and  composition  control.  Relaxing  the  structure  in  a Lennard- Jones  potential 

improved  the  isotropy,  the  density  uniformity  and  the  appearance  of  the  split 

31 

second  peak  in  the  radial  distribution  function.  For  a 20%  metalloid  alloy 
they  obtained  coordination  numbers  (see  Table  1)  and  a structure  factor  that 
compare  very  well  with  experiment.  For  a complete  evaluation,  their  partial 
radial  distribution  functions  should  be  compared  in  detail  with  the  experi- 
mental ones.  Considering  the  earlier  discussions  of  atomic  sizes,  this  pro- 
bably requires  a careful  choice  of  the  sphere  size  ratio  and  the  potential 
parameters. 

In  conclusion,  it  can  be  said  that  the  qualitative  aspects  of  the  Polk 
3 5 

model  ’ (metalloid  surrounded  by  metal  atoms  only,  short  metal- metalloid 

distance,  metalloid  coordination  number  smaller  than  that  of  the  metal)  are 

borne  out  by  the  recent  scattering  studies,  but  that  some  relaxation  of  the 

metal  skeleton  is  necessary.  The  idea  of  a metal  skeleton,  however,  besides 

being  quite  elegant,  remains  an  attractive  one  in  the  light  of  Turnbull's 
24 

scaling  relation.  It  would  therefore  be  interesting  to  check  whether  the 
hole  structure  formed  by  the  metal  atoms  in  the  binary  DRP's  can  be  related 
simply  to  that  of  the  single  atom  DRP. 
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type  of  TSRO  and  this  tends  to  keep  defects  localized.  In  an  amorphous 

structure,  this  requirement  is  relaxed  and  therefore  a variety  of  TSRO 

types  can  be  accommodated.  Because  of  this  flexibility,  it  seems  that 

localized  point,  line  and  planar  defects  are  not  likely  to  occur  as  such 

in  amorphous  systems,  but  probably  become  diffuse  by  breaking  up  in 

small  perturbations  of  the  TSRO  over  an  extended  region.  This  has 

43 

indeed  been  observed  for  artificially  created  "vacancies"  in  two-, 

44 

and  three-dimensional  amorphous  model  systems. 

4 7 

Based  on  positron- annihilation  studies  of  cold  rolled  metallic 
glasses,  the  same  can  probably  be  expected  for  a line  defect.  The 
analogue  of  a localized  crystalline  dislocation  created  instantaneously 
in  an  amorphous  model  system  would  probably  be  unstable.  Not  only 
would  its  core  become  very  diffuse,  but  the  displacements  and  climb 
resulting  from  its  interaction  with  local  stress  fields  and  partial 
vacancies  would  lead  to  a dislocation  line  that  is  contorted  on  an  atomic 
scale.  Since  under  the  influence  of  an  external  stress  these  segments 
will  move  in  an  uncorrelated  way,  the  resulting  deformation  process 
should  be  thought  of  as  one  governed  by  the  motion  of  individual  defect 
sites  rather  than  by  glide  of  a dislocation  line. 

The  most  fruitful  way  to  describe  structural  defects  in  amorphous 
systems,  therefore,  is  as  individual  sites  where  the  'ideal'  or  preferred 
local  SRO  is  perturbed.  In  general  they  include  both  deviations  from 
the  ideal  CSRO  (i.  e.  , wrong  nearest  neighbors)  or  ideal  TSRO  (e.  g.  , 
the  last  configuration  in  Fig.  5,  which  has  probably  a higher  energy 
than  the  first  one).  The  presence  of  non-ideal  SRO  facilitates  structural 
rearrangements  at  these  sites,  since  it  lowers  the  energy  difference 

r 
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between  the  initial  and  activated  state  of  the  rearrangement  process. 

The  various  atomic  transport  processes  (diffusion,  volume  relaxation, 
flow)  are  therefore  governed  by  the  concentration  and  mobility  of  these 
defects . 

43 

It  has  been  suggested  that  the  large  differences  that  have  been 
46 

observed  between  the  time  constants  of  some  of  the  atomic  transport 
processes  can  be  explained  by  distinguising  between  different  types  of 
defects  depending  on  the  type  of  rearrangement  they  allow  (nearest- 
neighbor  change  -»  diffusion;  long  range  stress  field  — volume  relaxation; 
local  shear  — flow). 
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IV.  STEADY  STATE  PLASTIC  FLOW 

1 . Survey  of  the  Experiments 

A convenient  way  of  summarizing  the  flow  data  for  a material  is 

48 

by  means  of  an  empirical  deformation  mechanism  map:  each  point 

on  the  map  represents  an  experiment  at  constant  stress  (j ) and 

temperature  (T),  and  contours  are  drawn  connecting  points  of  the  same 

steady  strain  rate  (y).  Figure  6 is  an  example  for  Pd-based  glasses. 

37-39  32  35 

Data  obtained  from  compression  and  creep  tests  ’ , or 

measurements  of  the  ultimate  tensile  strength  (the  three  highest  points 
of  Ref.  33)  are  probably  close  to  the  steady  state  values.  For  all  the 
other  tests  the  fracture  stress,  which  in  most  of  these  cases  coincides 
with  the  yield  stress,  has  been  plotted.  It  is  probably  too  low,  since 
steady  state  flow  has  been  forestalled  by  fracture.  If  more--and  more 
consistent- - data,  or  complete  constitutive  equations  were  available,  the 
resulting  deformation  mechanism  map  would  probably  resemble  the 
schematic  one  of  Fig.  7.  It  is  included  here  to  illustrate  the  basic 
modes  of  deformation:  homogeneous  and  inhomogeneous  flow. 

(a)  Homogeneous  Flow 

In  this  type  of  flow  each  volume  element  of  the  specimen  contri- 
butes to  the  strain.  In  a uniaxial  tensile  test,  fracture  occurs  at  large 
strains  and  after  extensive  necking.  This  flow  mechanism  operates  at 
low  stress  levels  at  all  temperatures  and  is,  except  near  the  boundary 
of  the  inhomogeneous  flow  region,  almost  always  a Newtonian  viscous 
(i.  e.  , y oc  t).  For  this  reason  the  shape  of  the  strain  rate  contours 


-22- 


reflects  that  of  the  viscosity- -temper ature  curve, 
arbitrarily  defined,  regions  can  be  distinguished: 


43 


Three,  somewhat 


(i)  T > (liquidus  temperature).  Only  one  set  of  viscometric 

measurements  on  a glass  forming  melt  (Au^^Si^  ^Ge^  4)  has  been 

reported.  50  The  viscosity  is  similar  to  that  of  other  high  temperature 

liquids:  q 10'  poise,  ('fluid')  and  changes  only  slowly  with 

temperature  as  indicated  by  the  horizontal  strain  rate  contours  in  Fig.  7. 

(ii)  T *a  T (glass  transition  temperature) . The  viscosity  is  by 
1 3 

definition  ~10  poise  in  this  region  and  hence  the  system  could  be 

called  "viscous".  A number  of  creep-type  viscosity  measurements 

Si  3 Z S Z 

have  been  reported  on  Au-based  and  Pd-based  ’ glasses.  Since 

in  this  temperature  range  the  rate  of  thermal  tr ansformations  seems 

to  scale  with  the  viscosity,  it  is  possible  to  evaluate  the  viscosity 

from  a set  of  calorimetric  ally  obtained  crystallization  and  glass 
52 

transition  data.  In  all  cases  the  data  could  be  represented  in  this 

temperature  range  by  a Fulcher- Vogel  type  equation:  r]  = A exp[B/(T- Tq)  ]. 

The  strain  rate  contours  on  Figs.  6 and  7 rise  very  steeply  with 

temperature,  i.  e.  , the  viscosity  has  a large  apparent  activation  energy. 

15 

(iii)  T < T . The  system  has  a viscosity  r\  > 10  poise  and  can 

therefore  be  called  "solid.  " Creep  measurements  in  this  region  are 

complicated  by  continuous  structural  relaxation  of  the  system  towards 

54  35 

higher  viscosities.  For  this  reason,  the  early  measurements 


on  Ni?5P25,  Co 


P and  Pdor.Sior>  should  be  treated  with  caution. 
75  Z5  oU  ZU 


For  example,  rerunning  the  measurements  on  Co-^P^  gave  a higher 
viscosity,  which  indicated  that  structural  changes  had  taken  place.  The 
low  activation  energy  observed  for  PdggSi^g  (0.  6 eV)  might  also 


r 

-23- 

increase  up  annealing,  as  indicated  by  the  diffusion  measurements  on 
46 

a similar  system.  In  order  to  perform  a truly  iso- configurational 

creep  test  it  is  necessary  to  anneal  the  specimen  for  a long  time  at  a 

temperature  higher  than  that  of  the  test,  as  demonstrated  by  Hadnagy 
55 

et  al.  in  their  load  relaxation  studies  of  a Fe^QNi^gPj^Bg  glass. 

They  found  that  the  plastic  flow  component  of  the  relaxation  had  a strain 

rate  sensitivity  exponent  (m  = d In  y/d  In  t)  of  about  4.  This  is 

35  54 

larger  than  what  had  been  measured  earlier  ’ (m  » 1),  but  similar 
to  Megusar  et  al.'s  recent  observations^  on  PdggSi^Q  (4<  m < 10). 

(b)  Inhomogeneous  Flow 

This  type  of  flow,  in  which  the  strain  is  localized  in  a few  very 

thin  shear  bands,  operates  at  high  stress  levels.  (See  Fig.  7)  In  this 

region,  the  strain  rate  contours  are  very  close  together,  i.  e.  , the 

strain  rate  is  very  stress  sensitive  (the  exponent  m is  large).  The 

flow  stress  is  almost  temperature  independent,  except  for  an  upturn  at 

61  62 

low  temperatures  especially  in  Fe-  and  Ni-based  glasses.  ’ Recent 

direct  observations  of  the  shear  bands  have  shown  that  they  are  ~50  A 
5 7 58 

thick  and  form  very  rapidly  (<  6.  7 msec.  I prior  to  fracture. 

Although  the  specimen's  total  plastic  strain  is  quite  small  (<  0.  2%), 
the  local  plastic  strain  in  the  shear  bands  is  very  large.  This  causes 
the  local  cross  section  of  the  specimen  to  decrease  until  the  stress 
concentration  is  large  enough  to  induce  fracture  along  the  shear  plane 
by  the  Taylor  instability  mechanism,  ^ which  produces  a vein-like 
pattern  on  the  fracture  surface. 

It  seems  clear  from  a number  of  considerations  that  this  inhomo- 


geneous flow  phenomenon  must  involve  a structural  change  which 
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produces  a local  softening  (i.  e.  , lowering  of  the  viscosity)  in  the  shear 
band:  (i)  such  a change  would  obviously  confine  all  shear  to  the  band; 

(ii)  the  shear  bands  occur,  depending  on  the  specimen  geometry,  at 
different  angles  with  the  tensile  axis  but  are  never  normal  to  it;  since 
fracture  occurs  along  the  shear  band  rather  than  normal  to  the  tensile 
axis  indicates  that  the  band  has  been  weakened;  (iii)  differential 
etching^’ ^ ancj  small  changes  in  the  diffraction  pattern*^  have  been 
observed. 

2.  Microscopic  Models 

(a)  Defect- Free  Flow:  The  Ideal  Shear  Strength 

The  ideal  shear  strength  (t.^)  is  the  stress  at  which  the 
structure  becomes  unstable:  it  shears  off  over  its  whole  cross-section 
at  once.  The  simple  Frenkel  calculation,  ^ which  assumes  a sinusoidal 
shape  for  the  potential  curve,  yields  t.^  = p/2ir  (p:  shear  modulus). 

It  is  probably  more  appropriate  for  amorphous  solids,  which  are 
mechanically  isotropic,  then  for  crystals,  where  subsidiary  minima 
in  the  potential  curve  corresponding  to  the  easiest  glide  directions  tend 
to  give  a lower  value  for  t.^. 

(b)  Defect  Controlled  Flow 

(i)  The  General  Flow  Equation.  As  discussed  in  the  previous 
section,  atomic  transport  processes  are  governed  by  the  concentration 
and  motion  of  the  appropriate  type  of  defect.  In  the  case  of  plastic  flow 
this  defect  is  a site  of  volume  Vq,  which  upon  rearrangement  produces 
a local  shear  strain  y^.  It  is  easy  to  show  that  the  macroscopic  strain 
rate  can,  quite  generally,  be  expressed  as: 
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i = nvo7oko 


where  n is  the  concentration  of  defect  sites  and  kg  the  frequency  of 
rearrangement  at  a site.  Chemical  rate  theory67’ 68  gives  an  expression 
for  kQ  of  the  form: 


vD  exp 


(-  w) sinh  (“W2) 


(21 


The  hyperbolic  sine  is  the  stress-dependent  driving  factor  and  AG'  is 
the  activation  free  energy  for  rearrangement,  which,  in  general, 
contains  chemical  and  strain  energy  contributions.  vD  is  the  frequency 

of  atomic  vibration  (~  Debye  frequency). 

Expressions  of  the  type  of  eq.  (1)  have  been  used  for  a long  time 
to  describe  flow68  and  are  general  enough  to  include  the  case  of  dis- 
location glide  in  crystals.  For  example,  in  a simple  cubic  lattice 
(lattice  parameter  a),  n = p/l  (p:  dislocation  density;  I : average 
dislocation  length),  vQ  = a2f,  yQ  = b/a  (b:  Burgers  vector),  kQ  = v/a 
(v:  dislocation  velocity).  Equation  ( 1)  then  becomes  y = pbv,  the 
Orowan  equation.  Based  on  the  macroscopic  analogy  between  the 
surface  steps  produced  by  the  shear  bands  in  amorphous  metals  and 
the  slip  bands  formed  during  plastic  deformation  of  crystals,  a number 

of  workers6^  7^  have  invoked  the  presence  of  localized  dislocation  lines 

. . 47  ....  . 

to  explain  inhomogeneous  flow.  Subsequent  investigations,  however, 
have  failed  to  produce  evidence  for  their  existence  and  it  will  be  shown 
later  on  that  the  formation  shear  bands  can  be  explained  as  a macro- 
scopic instability  effect,  without  recourse  to  dislocation  lines.  The 
present  approach,  based  on  local  shear  produced  by  individual  defect 
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sites,  is  more  consistent  with  the  structural  observations  and  provides 
a more  unified  view  of  the  deformation  process  since  it  applies  to  both 
homogeneous  and  inhomogeneous  flow. 

(ii)  Homogeneous  Flow,  T > ( 'Fluid') . An  adequate  character- 

ization of  the  viscosity  of  liquids  at  high  temperatures  has  been 

67  7 Z 

available  for  a long  time.  ’ The  mechanism  involves  switching  of 
individual  atoms  (Vq  atomic  volume  fi,  Yq  w 1).  The  number  of 
defect  sites  is  of  the  same  order  of  magnitude  as  the  number  of  atoms 
and  therefore  changes  little  with  temperature.  For  the  low  stress 
levels  of  fluid  flow:  sinh(T  YqVq/IcT)  sa  tYqVq/IcT,  i.  e.  : the  flow  is 
Newtonian.  The  activation  free  energy  AG'  is  very  small,  probably 
because  at  these  high  defect  concentrations  the  structure  is  flexible 
enough  to  make  the  elastic  energy  involved  in  switching  negligible.  For 
this  reason  the  temperature  dependence  of  the  viscosity  is  small. 

(iii)  Homogeneous  Flow.  TwTg  ('Viscous ').  In  this  region 
the  change  of  \ with  temperature  is  dominated  by  the  defect  concen- 
tration n.  The  SRO  in  the  system  increases  rapidly  with  decreasing 
temperature,  and  this  causes  a corresponding  decrease  of  the  defect 
concentr  ation. 

7-2 

In  the  free  volume  theory,  developed  by  Turnbull  and  Cohen 
for  monatomic  systems,  a defect  site  is  defined  as  a local  density 
fluctuation  which  creates  a hole  large  enough  to  allow  a neighboring 
atom  to  jump  into  it  (i.  e.  , v^  as  fl  , Yq  rs  1).  The  defect  concentration 
is  then 
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where  is  a geometrical  overlap  factor  of  order  unity,  v*  the 
minimum  hole  size  (~  ion  core  volume  for  metals)  and  v^  the  average 
free  volume  per  atom  (defined  as  the  difference  between  the  actual 
atomic  volume  (H)  and  that  of  a reference  system  with  ideal  SRO  (Oq); 
the  reference  system  for  hard  spheres  is  the  DRP).  At  constant 
pressure,  the  average  free  volume  can  be  expressed  as  v^  = afT-T^)!^, 
where  a is  an  average  coefficient  of  thermal  expansion  and  Tq  the 
equilibrium  temperature  of  the  reference  system.  Combined  with 
eq.  (3),  which  is  the  dominating  term  in  eq.  (1),  this  leads  to  a Fulcher- 
Vogel  type  expression  for  the  viscosity. 

Since  a is  not  a constant,  however,  and  measurements  of  it  are 
rarely  available,  theoretical  calculations  of  the  thermal  expansion  are 

-J  A 

useful.  This  has  been  done  by  Ramachandrarao  et  al.  , who  used  the 
hole  theory  of  liquids  to  derive  an  expression  for  v^  which  varies 
exponentially  with  inverse  temperature.  This  allowed  them  to  fit  the 
viscosity  data  for  a member  of  glassforming  melts  over  the  entire 
temperature  range  T.  - T something  which  cannot  be  done  with  a 

* fe 

single  Fulcher- Vogel  expression. 

An  alternative  approach  to  the  calculation  of  the  defect  concen- 

?C 

tration  is  the  configurational  entropy  model  of  Adam  and  Gibbs,  in 
which  a defect  site  is  defined  as  a subsystem  with  a configurational 

nL 

entropy  exceeding  some  critical  value  S*.  Chen  has  applied  this 
model  to  the  flow  of  glassy  metals.  It  appears  that  under  certain 
conditions  the  two  models  should  reduce  to  the  same  form.  The  con- 
figurational entropy  approach,  however,  is  more  phenomenological  in 
nature  and  does  not  seem  to  provide  as  much  insight  into  the  atomistic 
details  of  the  local  shear  mechanism  as  the  free  volume  model  does. 
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Since  the  free  volume  model  implies  Vq  "Yq  % 0 one  would 
expect  that  at  the  experimental  stress  levels  sinh(T  Vq  -y^/kT)  « TV^-Y^/kT 

and  that  the  flow  would  be  Newtonian.  Chen  and  Goldstein's  measure- 
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ments  on  Pdyy  j-Cu^Si^  glasses  confirm  this  for  high  viscosities 
(10*2  - lO1^  poise),  but  at  lower  viscosities  ( 10 1 1 - 10*2  poise)  they 
report  non- Newtonian  behavior,  corresponding  to  a value  of  v^  as 
high  as  100(1.  This  result  seems  somewhat  inconsistent  in  view  of  the 
Newtonian  behavior  at  both  lower  and  higher  viscosities,  and  it  remains 
as  yet  unexplained. 

(iv)  Homogeneous  Flow,  T < Tg  ( 'Solid').  In  order  to  perform 
meaningful  measurements  in  this  region  it  is  necessary  to  stabilize  the 
structure  by  annealing.  Since  such  a system  is  frozen  in  one  particular 
configuration,  one  might  at  first  expect  that  according  to  the  free  volume 
model  the  defect  concentration  would  remain  constant  at  the  level  defined 
by  the  fictive  temperature.  However,  it  should  be  realized  that  even  in  a 
frozen-in  configuration,  v^  still  decreases  with  decreasing  temperature 
because  of  the  non-configurational  (i.  e.  , uniform)  part  of  the  thermal 
expansion.  At  the  same  time  v*  increases  because  of  the  increasing 
effective  hard  core  diameter  of  the  metal  ion  (defined  as  the  distance 
at  which  the  repulsive  potential  energy  is  kT  above  the  minimum). 
Equation  (3)  shows  that  these  two  effects  result  in  a continued  decrease 
in  defect  concentration  with  decreasing  temperature.  Combined  with  a 
contribution  from  AG1,  which  in  this  model  is  mainly  chemical  in 
origin,  this  would  predict  a fairly  large  observed  activation  energy  for 
the  iso- configurational  viscosity.  Such  data  are  not  available  yet,  but 
the  large  activation  energy  observed  in  diffusion  measurements46  on 
annealed  specimens  points  in  the  same  direction. 
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Recently,  Argon  has  proposed  a somewhat  different  approach 
to  this  type  of  flow.  In  the  three  volume  description  the  defect  con- 
centration  is  small  (n  < 10  /G),  but  the  rearrangement  at  each  site 

is  relatively  easy  (AG1  contains  only  a chemical  contribution).  In 
Argon's  model,  the  defect  sites  are  larger  (v^  *»  50  G)  but  the  local 
shear  strain  is  smaller  (y^  rs  0.  1);  the  defect  concentration  is  larger 
(n  Vq  % 1)  but  the  rearrangements  are  more  difficult  (AG1  contains 
an  elastic  strain  energy  and  a shear  resistance  term).  The  qualitative 
results  from  the  two  models  are  therefore  likely  to  be  similar. 

The  models  differ  quantitatively  in  their  predictions  of  the  strain 
rate  sensitivity  exponent.  From  eq.  (1)  and  (2)  follows: 


m 


d In  y 
3 In  t 
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Observations  have  been  made  for  stresses  around  j = 5 x 10  Nm 
and  T = 500K.  Argon's  model  (yQ  vo  **  predicts  m « 5,  in 
accordance  with  the  observations.  The  free  volume  model  (Yqvo  ***  ^ 
predicts  m = 1,  and  needs  therefore  some  generalization  which  allows 
the  rearrangement  of  more  atoms  per  defect  site. 


(v)  Inhomogeneous  Flow.  As  discussed  above,  formation  of  the 
shear  bands  involves  a structural  change  that  leads  to  a lowering  of  the 
viscosity  in  the  band.  In  the  free  volume  model,  the  viscosity  is 
governed  by  the  defect  concentration,  and  hence  a decrease  in  viscosity 
must  correspond  to  an  increase  in  the  average  free  volume  v^. 

Following  a general  suggestion  by  Polk  and  Turnbull,  50  Spaepen49  has 
shown  that  at  high  stress  levels  the  free  volume  in  the  band  is  increased 
as  a result  of  a dynamic  equilibrium  between  stress-driven  creation  of 
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free  volume  (atoms  pushing  their  neighbors  aside)  and  its  annihilation 
by  diffusional  rearrangements  (the  system  tries  to  relax  back  to  its 
original  state).  The  resulting  expression  for  the  defect  concentration 
does  not  depend  on  the  original  structure,  but  is  a strong  function  of 
the  stress  and  the  diffusion  rate: 


In  n 


S v* 
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where  S = 2p(l + v)/3(l  - v)  and  n^  the  number  of  jumps  necessary  to 
annihilate  an  amount  of  free  volume  equal  to  v*  (l<n^<  10). 

The  model  can  predict  the  boundary  line  between  the  homogeneous 
and  inhomogeneous  flow  regions  (see  Fig.  8)  in  terms  of  reasonable 
physical  parameters.  Since  the  defect  concentration  n is  the  dominant 
term  in  eq.  (1)  and  a strong  function  of  the  stress,  the  model  predicts  a 
very  large  strain  rate  sensitivity  coefficient  m,  in  accordance  with 
observations.  An  estimate  of  m can  easily  be  made  by  differentiating 
eq.  (5)  as  a function  of  In  t . Insertion  of  typical  values  (S  = 5.  12x  10*9 
Nm'2,  T = 6xl08Nm"2,  Q = 1.4x  10'29m3,  v*/nD=0.26fi,  T = 300K) 

yields  m = 85. 

79 

Argon  has  recently  applied  this  idea  to  a model  involving 
"small  dislocation  loop-like"  defect  sites  (vQ  20  0,  Vq  = 1).  In  his 
model  the  free  volume  affects  the  viscosity  by  lowering  the  shear 
resistance  term  in  AG1.  It  predicts  the  homogeneous-inhomogeneous 
transition,  but  the  idea  of  a "small  dislocation  loop"  should  not,  for 
topological  reasons,  be  taken  to  literally.  More  importantly,  however, 
his  analysis  shows  rigorously  that  localization  of  the  shear  in  bands  is 
a direct  result  of  the  softening  caused  by  the  dynamic  excess  of  free  volume  . 
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